Parametric characterisation of a chaotic attractor using two scale Cantor measure 
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A chaotic attractor is usually characterised by its multifractal spectrum which gives a geometric 
' measure of its complexity. Here we present a characterisation using a minimal set of independant 

parameters which are uniquely determined by the underlying process that generates the attractor. 
The method maps the f{a) spectrum of a chaotic attractor onto that of a general two scale Cantor 
measure. We show that the mapping can be done for a large number of chaotic systems. In order 
to implement this procedure, we also propose a generalisation of the standard equations for two 
scale Cantor set in one dimension to that in higher dimensions. Another interesting result we have 
obtained both theoretically and numerically is that, the f(a) characterisation gives information only 
upto two scales, even when the underlying process generating the multifractal involves more than 
, two scales. 

. PACS numbers: 05.45.Ac, 05.45.Tp, 05.45.Df 
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I. INTRODUCTION 

The existence of a multifractal measure for any system most often indicates an underlying process generating it, be 
it multiplicative or dynamic. In the context of chaotic attractors arising from dynamical systems, their multifractal 
measure result from a time ordered process, which maybe an iterative scheme or a continuous flow |l[ . The description 
. of the invariant measures in terms of D q [2| or /(a) [3|, however, provides only a characterisation of their geometric 
complexity. Feigenbaum et.al 0, Q and Amritkar and Gupte p have shown that it is also possible to get the 
dynamical information in some specific cases by inverting the information contained in a multifractal measure using 
a thermodynamic formalism. 

In this paper, we seek to get a characterisation of a chaotic attractor in terms of the underlying process that 
generates it. It appears that the process of generation of a multifractal chaotic attractor is similar to that of a typical 
Cantor set (where measure reduces after each step), with the dissipation in the system playing a major role. We show 
this specifically below using the example of Cat map which is area preserving. But a key difference is that, for chaotic 
attractors, the nature of this reduction is governed by the dynamics of the system. This implies that if the D q and 
f(a) curves of a chaotic attractor are mapped onto that of a model multiplicative process, one can derive information 
about the underlying process that generates the strange attractor, provided the mapping is correct. Here we try to 
implement this idea using an algorithmic scheme and show that this gives a set of parameters that can be used to 
characterise a given attractor. 

A similar idea to extract the underlying multiplicative process from a multifractal has been applied earlier by 
Chhabra et.al Q- In order to make this inversion process successful, one needs to take into account two aspects, 
namely, the type of process Q (whether L process, P process or LP process) and the number of scales involved 
(whether two scale or multi scale). Chhabra et.al 0] have shown that different multiplicative processes with only three 
independant parameters produce good fits to many of the observed D q curves. Thus the extraction of underlying 
multiplicative process, based solely on the information of D q curve, is nonunique and additional thermodynamic 
information is needed for the inversion process. 
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But the problem that we address here is slightly different. In our case, the model multiplicative process is fixed as 
a general two scale Cantor set which is the simplest nontrivial process giving rise to a multifractal measure. We then 
scan the whole set of parameters possible for this process (which include the L process, P process and LP process) 
and choose the statistically best fit D q curve to the D q spectrum computed for the attractor from the time series, 
which is then used to compute the final f(a) spectrum. In this way, the f{a) spectrum of a chaotic attractor gets 
mapped onto that of a general two scale Cantor set. We show that the mapping can be done for a large number of 
standard chaotic attractors. The resulting parameters can be considered to be unique to the underlying process that 
generates the attractor, upto an ambiguity regarding the number of scales involved. 

The success of this procedure also implies that the D q and f(a) spectrum of a multiplicative process involving 
more than two scales also can be mapped onto that of a two scale Cantor set. We prove this theoretically as well as 
numerically in Sec. IV, by taking Cantor sets with more than two scales. This, in turn, suggests that though the f(a) 
spectrum has contributions from all the scales involved in the generation of a multifractal, the information contained 
in an f(a) spectrum is limited only upto two scales. In other words, given an /(a) spectrum, one can retrieve only the 
equivalent two scales which are different from the actual scales. Thus, while Chhabra et.al 0] argues that additional 
information is needed to extract the underlying multiplicative process, our result indicate that the /(a) formalism 
itself is unable to extract more than two scales. 

The motivation for using a Cantor set to characterise the multifractal structure of a chaotic attractor comes from 
the fact that some well known chaotic attractors are believed to have underlying Cantor set structure. For example, 
it has been shown [8( that in the x — y plane corresponding to z = (r — 1) of the Lorenz attractor, a transverse cut 
gives a multi fractal with Cantor set structure. Even the chaotic attractor resulting from the experimental Raylcigh- 
Bernard convection holds a support whose transverse structure is a Cantor set @. These Cantor sets are known to 
be characteristic of the underlying dynamics that generate the attractor. 

A more general arguement to support the above statement is by using the concept of Kolmogorov entropy. Kol- 
mogorov entropy can be obtained by a successively fine partition of the attractor in a hierarchical fashion. Going from 
one partition to the next gives one set of scales as shown in @ . These can be treated as scales of higher dimensional 
Cantor sets. In general, there can be several scales. But the f(a) curve appears to be determined by only two scales. 

In order to implement our idea, the first step is to compute the D a spectrum of the chaotic attractor from its time 
series. This is done by the standard delay embedding technique [1 01 ] . but by extending the nonsubjective scheme 
recently proposed by us [Tlj for computing D^. The D q spectrum is then fitted by a smooth D q curve obtained 
by inverse Legendre transformation equations [l2l Il3j of the /(a) curve for a general two scale Cantor set. The 
statistically best fit curve is chosen by changing the parameters of the fit from which, the /(a) curve for the time 
series is evaluated along with a set of independant parameters characteristic of the Cantor set. This procedure also 
gives a couple of other interesting results. For example, we are able to propose a generalisation of the standard 
equations of two scale Cantor set for higher dimensions. Moreover, we explicitely derive the equations for D q and 
f(a) spectrum of a three scale Cantor set and show that they can be exactly mapped onto that of a two scale Cantor 
set. 

Our paper is organised as follows: The details of our computational scheme are presented in Sec. II and it is tested 
using time series from the logistic map and different Cantor sets with known parameters in Sec. III. In Sec. IV, the 
f(a) spectrum of Cantor sets with more than two scales is considered both theoretically and numerically. Sec.V is 
concerned with the application of the scheme to standard chaotic attractors in higher dimensions. The conclusions 
are drawn in Sec. VI. 

II. NUMERICAL SCHEME 

As the first step, the spectrum of generalised dimensions D q are computed from the time series using the delay 
embedding technique [lcj . For a given embedding dimension M, the D q spectrum are given by the standard equation 

D q , — hm ^ °^ (1) 
q q-1 log R W 

where C q (R) represents the generalised correlation sum. In practical considerations, D q is computed by taking the 
slope of log C q (R) versus logi? over a scaling region. In our scheme, the scaling region is computed algorithmically 
[111 ] for each D q using conditions for R m in and R ma x and the spectrum of D q for q in the range [—20, 20] is evaluated 
with an error bar. 

Assuming that the corresponding f(a) curve is a smooth convex function, we seek to represent it using the standard 
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FIG. 1: The upper panel shows the D q values (points) and its best fit curve of the strange attractor at the period doubling 
accumulation point of the logistic map, computed from a time series with 30000 data points. The lower panel shows the f(a) 
spectrum computed from the best fit curve (dashed line) along with the theoretically predicted curve (continuous line). The 
agreement between the two is evident. 



equations [HdH of a and f(a) for the general two scale Cantor set 

rlogpi + (1 - r)logp 2 

O! = 

rlog/i + (1 — r) log/ 2 



(2) 



= rlogr + (1 -r) log(l - r) 
1 rlog/ 1 + (l-r)logZ 2 (> 

where l\ and li are the rescaling parameters and p\ and pi are the probability measures with p2 = (1 — pi). Thus there 
are three independent parameters which are characteristic of the multiplicative process generating a given f{a) curve. 
Here r is a variable in the range [0, 1], with r — > corresponding to one extreme of scaling and r — > 1 corresponding 
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to the other extreme. Taking logp2/log/2 > logpi/logli, as r — > 0, we get 

logl 2 

and as r — > 1 

= logPi 
log/i 



(4) 



(5) 



By inverting Eqs. © and ([3]) and using the standard Legendre transformation equations [l2l [l3j connecting a and 
/(a) with q and Z? q , we get 

q = /(«) (6) 

da 



m - /(« ) 

(3-1) 



Changing the variable rj = 1/r, ^ and reduce to 

log Pi + (v~ l)logP2 

OL = 

logli + (77 - l)logl 2 



(8) 



, = (y- i)iog(>7- 1) -yiogy (q) 

1 log/ 1 + (77-1) log Z 2 K) 



Differentiating ([5]) and ^ with respect to r\ and combining 

d£_ = (logZi(log(T?- 1) -log?7) + log?2 log??) 
da (logZilogp 2 -logZ 2 logpi) 

Using eq. ^ and changing back to variable r 

df log h log(l - r) - log l 2 log r 



da logZilog(l-pi)-logZ 2 logpi 



(10) 



(11) 



Eqs. pip and ([7]) give both g and D q as functions of the three independent parameters li, Z 2 and pi. 

For a given set of parameters, the D q curve is determined by varying r in the range [0, 1] and fitted with the 
computed D q values from the time series. The procedure is repeated by changing the values of p\ in the range [0, 1] 
and for each pi, scanning the values of l\ and I2 with the condition that both < 1. A statistical x 2 fitting is 
undertaken and the best fit curve given by the % 2 minimum is chosen. The complete f(a) curve is derived from it 
along with the complete set of parameters pi, Zi, la, a m in and a m ax, for a particular time series. 



III. TESTING THE SCHEME 



In order to illustrate our scheme, we first apply it on standard multifractals where the f(a) curve and the associated 
parameters are known exactly. In all the examples discussed in this paper, 30000 data points are used for the analysis. 
The first one is the time series from the logistic map at the period doubling accumulation point. The D q spectrum 
is first computed using Eq. |T]) (with M = 1), for q values in the range [—20, +20]. The computation is done taking 
a step width of Aq — 0.1. Choosing pi = 0.5,a m i n = D20 and a max — -D-20 as input parameters, the D q curve is 
computed from the above set of equations and fitted with the D q values. The procedure is repeated by scanning pi 
in the range [0, 1] in steps of 0.01. For each pi, a m i n and a max (which in turn determine li and l 2 ) are also varied 
independantly over a small range. The best fit D q curve is chosen as indicated by the x 2 minimum. Since the error 
in D q generally bulges as q — > —20, the error bar is also taken care of in the fitting process. The D q values computed 
from the time series and its best fit curve are shown in Fig. [TJ 

The complete f(a) spectrum for the time series is computed from the best fit D q curve. To make a comparison, the 
spectrum is also determined from Eqs. and ^ using the known values of p\, l\ and Z 2 for the logistic map, namely, 
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FIG. 2: The D q values (dashed lines) of four different two scale Cantor sets computed from their time series along with the 
numerical fit (continuous line) in each case. 



Pi = 0.5, l\ = 0.158(l/a|.) and I2 — 0.404(1/0?) where qj? is Feigenbaum's universal number. Both the curves are 
also shown in Fig. [1] The three parameters derived using our scheme are p\ = 0.5, Zi = 0.146 and I2 = 0.416 which 
are reasonably accurate considering the finiteness of the data set. 

As the second example, we generate time series from four Cantor sets using four different sets of parameters as 
given in Table |H Fig. [5] shows the computed D q values along with the best fit curves in all the four cases. Note that 
the fit is extremely accurate for the whole range of q in all cases. The corresponding f{a) curves, both theoretical and 
computed from scheme are shown in Fig. [3J The parameter values derived from our scheme in the four cases are also 
given in Table [I] for comparison. It is clear that the scheme recovers the complete f(a) spectrum and the parameters 
reasonably well. In order to convince ourselves that the scheme does not produce any spurious effects, we have also 
applied it to a time series from a pure white noise. The D q versus q curve for white noise should be a straight line 
parallel to the q axis with Dq — M. The corresponding f(a) spectrum would be a 5 function which has been verified 
numerically. 

From the numerical computations of two scale Cantor sets, we also find the following results: While the end points 
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FIG. 3: The f{a) curves for the Cantor sets shown in Fig. [2] computed from the best fit D q curve (dashed line) along with 
the theoretical curves (continuous line). The parameter values used for constructing the Cantor sets agree very well with those 
derived using our scheme in all the cases (see Table Q}. 



of the spectrum, a m in and a m ax, are determined by the ratios logpi/logZi and logp2/log/2, the peak value Do is 
determined by only the rescaling parameters l\ and As (Zi + 1%) increases (that is, as the gap length decreases), D 
also increases and D — > 1 as (li + 1 2 ) — > 1. In this sense, the gap length also influences the f{a) spectrum indirectly. 
We will show below that this is not true in the case of three scale Cantor set where we miss some information regarding 
the scales. We also find that as the difference between a rn i n and a max increases (that is, as the spectrum widens), 
more number of data points are required, in general, to get good agreement between theoretical and numerical f(a) 
curves. 
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TABLE I: Comparison of the parameters used for the generation of the different Cantor sets discussed in the text with those 
computed by applying our numerical scheme. Close to 30000 points are used for computation in all cases. 



Cantor set No. 



Parameters used 



Parameters computed 



Cantor set 1 


Pi 


= 0.60, h = 


0.22, h = 


0.48 


Pi 


= 0.58, 


h -■ 


= 0.21, 


h ■ 


= 0.49 


Cantor set 2 


Pi 


n A o 7 

= U.42, Li = 


U.22, 12 — 


U.OI 


Pi 


= U.4o, 


7 

tl : 


n o a 
= U.24, 


h - 


= U.D( 


Cantor set 3 


Pi 


= 0.66, h = 


0.18, l 2 = 


0.62 


Pi 


= 0.69, 


h -■ 


= 0.19, 


h ■ 


= 0.64 


Cantor set 4 


Pi 


= 0.72, h = 


0.44, h = 


0.48 


Pi 


= 0.66, 


h 


= 0.39, 


h ■ 


= 0.52 


3Scale Cantor set 


Pi 


= 0.25, p 2 = 


-- 0.35, p 3 -- 


= 0.4 
















k 


= 0.12, l 2 = 


0.35, h = 


0.18 


Pi 


= 0.50, 


h ■ 


= 0.26, 


h ■ 


= 0.52 


4Scale Cantor set 


pi = 0.34, p 2 = 0.38, 


P3 = 0.16 


, p 4 = 0.12 
















h = 0.12, h = 0.25, 


Is = 0.18, 


U = 0.08 


Pi 


= 0.58, 


h 


= 0.30, 


h ■ 


= 0.57 



IV. MULTI SCALE CANTOR SETS 



In this section, we consider the /(a) spectrum of a Cantor set with more than two scales. First we show the 
numerical results using our scheme. For this, we first generate the time series for a general 3 scale Cantor set and 
compute its D q spectrum. The geometrical construction of a general 3 scale Cantor set is shown in Fig. 0J At 
every stage, an interval gets subdivided into three so that the set involves 3 rescaling parameters li, h, h and 3 
probability measures p\, P2, P3 as shown. The numerically computed D q spectrum for a typical 3 scale Cantor set 
(with parameters given in Table [T| is shown in Fig. [5] (upper left panel). The D q curve can be very well fitted by 
a 2 scale Cantor set and the complete f{a) spectrum for the 3 scale Cantor set is evaluated (lower left panel). We 
have repeated our computations for a 4 scale Cantor set as well and the results are also shown in Fig. [5] (right panel). 
In both cases, the parameters used for the construction of the Cantor sets and those computed by our scheme are 
given in Table [U Thus it is clear that the f(a) spectrum cannot pick up the full information about the various scales 
and probability measures. No matter how many scales are involved in the generation of the multifractal, the f(a) 
spectrum can be reproduced by an equivalent 2 scale Cantor set. 

We now derive explicite expressions for a and f{a) for a 3 scale Cantor set. We follow the arguements given in 
Halsey et.al Q, Sec.II-C-4. For the 3 scale Cantor set, one can write 



'l l 2 '3 



Expanding 



T(q,T,n)= 



m\\m2^{in — mi — 7712)! 



miq m 2 q (m-m 1 -m 2 )q 1 -m 1 T ,-m 2 T ]-(m-m 1 -m 2 )T _ -. 
Pi Pi P3 L l L 2 L 3 ~ 1 



mi,7n 2 

In the limit n — > 00, the largest term contributes. Hence we have 

dT 



dmi 

dT 
dm 2 



= 



Using the Stirling approximation and simplifying the above two conditions we get 

-logr + log(l - r - s) + q\og(p 1 /p 3 ) -rlog(h/l 3 ) = 



(12) 



(13) 



(14) 



(15) 



(16) 



- logs + log(l - r-s) +q\og(p 2 /p 3 ) -T\og(l 2 /k) = 



(17) 
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k h h 



Pi" PjP 2 P1P3 PlP 2 Pa' P 2 P 3 PlP 3 P 2 P 3 P/ 



1 H 11 11 11 1 2 11 11 11 1 2 

M 1 l i 2 i l 1 3 1 1 1 B i 2 1 2 1 3 i l 1 3 1 2 i 3 i 3 



FIG. 4: The construction of a general 3 scale Cantor set. 

where r = mi/n and s — mijn are free parameters. Also from Eq. (|13p . using a similar procedure, one can show that 

r logr— slogs -(1-r-s) log(l-r-s) + <7(rlogpi + s log ^2 + (1-r-s) logp 3 ) --r(r logii+slogi2 + (l— r— s) log/3) = 

(18) 

Combining Eqs. (fl"6f . (jTTj) and (p~8|) and eliminating r we get the following relations for g 

^ = log(W log((l - r - s)/r) - log(Zi/Z 3 ) log((l - r - s)/s) 
log(/i// 3 ) log(p 2 /p 3 ) - \og{l 2 /h) log(pi/p 3 ) 

__ log(Zi/Z 3 )(— r logr - s log s - (1 - r - s) log(l — r — s)) — (rlogZi + slog/2 + (1 — r — s) log/ 3 ) log((l — r — s)/r) 
(rlogZi + slogZ 2 + (1 - r - s) log/ 3 ) log(pi/p 3 ) - log(/i// 3 )(r logpi + slogp 2 + (1 - r - s)logp 3 ) 

These two equations for q can be used to obtain a relation between r and s. 

To compute the D q spectrum, vary r from to 1. For every value of r, the value of s that satisfies the Eqs. I|19p 
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FIG. 5: The D q and f(a) curves of a 3 scale and a 4 scale Cantor set computed using our scheme. The D q curves computed 
from the time series (dashed lines) are fitted with that of equivalent two scale Cantor sets (continuous lines) in both cases. The 
parameters used and computed by the scheme are compared in Table U 



and (|20p simultaneously is computed numerically, with the condition that < s < (1 — r). For every value of r and 
s, q and r can be determined using Eqs. (|19|) and (|16p . which in turn gives D q = r/(q — 1). 
The singularity exponent a is determined by the condition 



„TOi_TO2_(m-m 1 -m2) _ /"roiimj/tm-mi-mj) 
Pi Pi P3 ~ l'l l 2 l 3 



This gives the expression for a 



r logpi + s logp 2 + (1 - r - s) logp 3 
r log h + s log h + (1 — r — s) log ^ 3 

Similarly, the density exponent /(a) is determined by 

,/(«) 



n! 1 1 1 l 2 tg ) = mi\rri2'-[rn — m\ — m2)\ 



(21) 



(22) 



(23) 
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FIG. 6: The left panel shows the theoretical D q curve of a 3 scale Cantor set with parameters shown in the figure. The 
right panel shows the corresponding f(a) curve (continuous line) calculated theoretically using Eqs. (|22|l and (|24p . This is 
superimposed on the f(a) spectrum (points) computed using our numerical code by fitting 2 scale Cantor set. Close to 30000 
data points are used for the computations. 



which gives the following expression for f(a) 

= r log r + s log s + (1 - r - s) log(l - r - s) 
1{a> rlogZi + slog h + (l-r-s) log l 3 1 ' 

By varying r from to 1, the f{a) spectrum for a given 3 scale Cantor set can be determined theoretically. In Fig. [6l 
the theoretically computed D q and f(a) spectrum for a typical 3 scale Cantor set is shown. Along with the theoretical 
f(a) curve, we also show the numerical one (points) for the same Cantor set, computed using our scheme. Thus it is 
evident that the f(a) spectrum of a 3 scale Cantor set can be mapped onto that of a 2 scale Cantor set. Also, our 
numerical results on 4 scale Cantor set (Fig. O suggests that this mapping onto 2 scale Cantor set can possibly be 
extended for the f(a) spectrum of four or more scale Cantor sets. 



11 



1.5 



0.5 



0.5 



Henon Attractor 
(M = 1) 




Lorenz Attractor 
(M = 1) 




0.5 



1.5 



2.5 



a 



FIG. 7: The multifractal spectrum of the one dimensional projection of Henon and Lorenz attractors, which show the underlying 
Cantor set structure. Note that the peak values of the spectrum in both cases are equal to 1. 



CHARACTERISATION OF STRANGE ATTRACTORS 



Evaluating the ,f{a) spectrum of one dimensional sets is straightforward. But computing the spectra of even 
synthetic higher dimensional attractors is a challenging task. Generally, the f(a) spectrum for higher dimensional 
chaotic attractors is calculated taking only one dimension [ll| [l7| , which characterise the transverse self similar 
structure on the attractor equivalent to a Cantor set. In the resulting /(a) spectrum, the peak value (that is, D ) 
will be equal to 1, as the higher dimensional attractor is projected into one dimension. This is shown in Fig. [7] for 
Henon and Lorenz attractors and the results are consistent with the earlier results. 

In order to extend our scheme to higher dimensional strange attractors, their /(a) spectra are to be considered 
analogous to a two scale Cantor measure in higher dimensions. While the /(a) curve can be recovered using the correct 
embedding dimension M, the meaning of the parameters have to be interpreted properly. For a one dimensional Cantor 
set, pi is a probability measure while l\ and I2 are fractional lengths at each stage. Extending this analogy to two and 
three dimensions, p\ can still be interpreted as a probability measure for the two higher dimensional scales, say T\ 
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FIG. 8: The upper panel gives the D q spectrum, with error bar, of the Rossler attractor computed from the time series. 
The lower panel shows the accuracy of fitting with again the D q values (dashed line) and its best fit curve (continuous line) 
computed using the scheme. 



and Ti . These can be considered as fractional measures corresponding to area or volume depending on the embedding 
dimension M . In other words, p\ is a measure representing the underlying dynamics, while t% and r 2 correspond to 
geometric scaling. This gives an alternate description of the formation of a strange attractor if it is correlated to a 
higher dimensional analogue of the Cantor set. 

As discussed in Sec. II, for the one dimensional Cantor set, a m in and a max are given by Eqs. ^ and |5]) in terms of 
the parameters. For p\ = p2 and l\ = I2, otmin — &max < 1 and the set becomes a simple fractal with a = /(a) = Dq, 
the fractal dimension. Extending this analogy to higher dimensions, we propose that Eqs. (0]) and © are to be 
modified as 

a ma x = M- (25) 

logr 2 
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FIG. 9: The f(a) curve of the Rossler attractor computed from the best fit D q curve shown in Fig. [8] 



and 



Otmin — M 



log Pi 

logn 



(26) 



As in the case of one dimensional Cantor sets, for p\ = p2 and t± = t%, a max = ot-min < M and the set is again a 
simple fractal with fractal dimension D = a = f(a). Rewriting the above equations, 



logp 2 logp 2 



and 



logfaVJ") log/ 2 

log Pi _ log Pi 
log(r 1 1 /M) " i ogh 



(27) 



(28) 
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TABLE II: The complete set of parameters computed using our scheme for six standard chaotic attractors. See [l8l ]. for example, 
for the details of the chaotic systems and the values of the parameters used. 



Attractor 




Qmax 


D 


Pi 




T2 


Rossler attractor 
(a = & = 0.2,c = 7.8) 


1.46 ±0.02 


3.39 ±0.14 


2.31 ±0.02 


0.65 


0.42 


0.41 


Lorenz attractor 
(cr = 10, r = 28, b = 8/3) 


1.38 ±0.03 


3.71 ±0.12 


2.16 ± 0.04 


0.50 


0.22 


0.57 


Ueda attractor 
(k = 0.05, A = 7.5) 


1.73 ±0.05 


3.78 ±0.13 


2.62 ± 0.06 


0.64 


0.46 


0.44 


Duffing attractor 
(b = 0.25, A = 0.4, Q, = 1) 


1.84 ± 0.04 


3.59 ± 0.08 


2.78 ± 0.04 


0.81 


0.71 


0.25 


Henon attractor 
(a = 1.4,6= 0.3) 


0.96 ±0.02 


2.27 ±0.08 


1.43 ±0.03 


0.50 


0.24 


0.54 


Tinkerbell attractor 
(a = 0.9, b = -0.6, c = 2, d = 0.5) 


0.83 ± 0.02 


3.43 ±0.12 


1.65 ±0.03 


0.60 


0.29 


0.58 



Replacing l\ and 1% by r\^ M and t^ M in Eqs. ([2]) and ([3]), the defining equations for the two scale Cantor set in M 
dimension can be generalised as 

a = M[r logpi + (1 -r) \ogp 2 ] ^ 
rlogri ± (1 - r) logr 2 



= M[rlogr±(l-r) log(l - r)] 
' r logTi ± (1 - r) logr 2 

Just like l\ + 12 < 1 for one dimensional Cantor set, we expect T\ ± r 2 < 1 in M dimensions. This is because, the 
measure keeps on reducing after each time step due to dissipation and t\ and t 2 represent the fractional reduction in 
the measure for the two scales. It should be noted that since, in general, different scales apply in different directions, 
T\ and r 2 should be treated as some effective scales in higher dimension. 

We now check these results using the time series from a standard chaotic attractor, namely the Rossler attractor for 
parameter values a = b = 0.2 and c = 7.8 with 30000 data points. Fig. [5] shows the D q spectrum computed from the 
time series taking M = 3, along with the best fit curve applying our scheme. The fit is found to be very good for the 
whole range of q values. The complete f{a) spectrum computed from the best fit D q curve is shown in Fig. [51 The 
scheme also calculates the three parameters as p\ = 0.65, t% = 0.42 and t 2 = 0.41 so that t% + t 2 = 0.83 < 1. Thus 
one can say that if the /(a) spectrum of the Rossler attractor is made equivalent to a two scale Cantor set in three 
dimension, the resulting probability measures are 0.65 and 0.35 and rescaling parameters 0.42 and 0.41. Interestingly, 
it appears that the Rossler attractor is generated by a P process rather than a LP process. 

The scheme has also been applied to several standard chaotic attractors in two and three dimensions. The /(a) 
spectrum are shown in Fig. [TU] for four of them, while the complete set of parameters for six standard chaotic 
attractors are given in Table [TXJ The error bars given for a max , a m i n and Dq are those reflected from the computed 
D q values. In a way, the two sets of parameters given above, that is pi,T%, r 2 and a max , a m i n , Dq, can be considered 
as complementary to each other. While the former contain the finger prints of the underlying process that generate 
the strange attractor (the extent of stretching and folding and redistribution of measures at each time step), the 
latter characterises the geometric complexity of the attractor once it is formed. Both can be independantly used to 
differentiate between chaotic attractors formed from different systems or from the same system for different parameter 
values. The former may be more relevant in the case of chaotic attractors obtained from experimental systems. 

Finally, we wish to point out that dissipation is a key factor leading to the multifractal nature of a chaotic attractor. 
To show this, we consider a counter example, namely, that of Cat map which is area preserving. The fixed points 
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FIG. 10: The multifractal spectrum of four standard chaotic attractors computed by applying our algorithmic scheme. The 
parameter values used for generating the time series are also shown. 



of the Cat map are hyperbolic, which are neither attractors nor repellers and the trajectories uniformly fill the phase 
space as time t — > oo. Its D q spectrum computed from the time series is found to be a straight line as shown in 
Fig. [Til just like that of a white noise. The corresponding f(a) curve is a 6 function, also shown in Fig. [TTJ Since 
ctmin = ctmax, a two scale fit gives the parameters as p\ = 0.5, T\ = 0.49 and T2 = 0.51. Thus the Cat map attractor 
turns out to be a simple fractal rather than a multifractal. 



VI. DISCUSSION AND CONCLUSION 



In this paper, we show that a chaotic attractor can be characterised using a set of three independant parameters 
which are specific to the underlying process generating it. The method relies on a scheme that maps the f(a) spectrum 
of a chaotic attractor onto that of a general two scale Cantor set. The scheme is first tested using one dimensional 
chaotic attractors and Cantor sets whose /(cv) curves and parameters are known and subsequently applied to higher 
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FIG. 11: The D q curve of the Cat map (left panel) computed using our scheme which is a straight line parallel to the q axis just 
like that of a white noise. The corresponding f(a) curve is shown in the right panel, which is 8 function with a = /(a) = 2. 



dimensional cases. 

In the scheme, the D q spectrum of a chaotic attractor is compared with the D q curve computed from a model 
multiplicative process. Similar idea has also been used to deduce certain statistical characteristics of a system and 
infer features of the dynamical processes leading to the observed macroscopic parameters. One such example has been 
provided earlier by Meneveau and Sreenivasan [191 ] in the study of energy dissipation rate in fully developed turbulent 
flows. By comparing the experimental D q data with that of a two scale Cantor measure, they have shown that the 
dynamics leading to the observed multifractal distributions of the energy dissipation rate can be well approximated 
by a single multi step process involving unequal energy distribution in the ratio 7/3. 

Usually, a multifractal is characterised only by the range of scaling involved [a m i n ,a max ], which roughly repre- 
sents the inhomogeinity of the attractor. So the set of parameters computed here seems to give alternative way of 
characterising them. But we wish to emphasize that the information contained in these parameters is more subtle. 
For example, once these prameters are known, a m in and a max can be determined using Eqs. (|27[) and (|28[) . Thus 
by evaluating pi, l\ and I2, we get additional information regarding the dynamic process leading to the generation 



17 



of the strange attractor. Moreover, these parameters can also give indication as to why the degree of inhomogeinity 
varies between different chaotic attractors. As is well known from the srudy of Cantor sets, the primary reason for the 
increased inhomogeinity is the wide difference between the rescaling measures l\ and ?2- Looking at the parameter 
values, rescaling measures t\ and t-i are very close for Rossler and Ueda attractors which appear less inhomogeneous, 
while that for Lorenz and Duffing are widely different making them more inhomogeneous with two clear scrolls. 

Another novel aspect of the scheme worth commenting is the use of two scale Cantor measures in higher dimension 
as analogues of chaotic attractors. Eventhough such objects are not much discussed in the literature, one can envisage 
them, for example, as generalisation of the well known Sierpinsky carpets in two dimension or the Menger sponge 
in three dimension. But a key difference between these and a chaotic attractor is that the rescaled measures are 
not regular in the generation of the latter. Recently, Perfect et. al [2(| present a general theoretical framework for 
generating geometrical multifractal Sierpinsky carpets using a generator with variable mass fractions determined by 
the truncated binomial probability distribution and to compute their generalised dimensions. It turns out that the 
chaotic attractors are more similar to multifractals generated in higher dimensional support, such as, fractal growth 
patterns and since the rescaled measures are irregular, a one dimensional measure such as l\ = t\^ M need not have 
any physical significance. 

Finally, for a complex chaotic attractor in general, the redistribution of the measures as it evolves in time can take 
place in more than two scales. Thus it appears that a characterisation based on only two scales is rather approximate 
as we tend to lose some information regarding the other scales involved. But we have found that the D q and /(a) 
curves of a multi scale Cantor set can be mapped onto that of an equivalent two scale Cantor set. These two scales 
may be functions of the actual scales involved and may contain the missing information in an implicit way. Thus, 
an important outcome of the present analysis is the realisation that the dynamical information that can be retrieved 
from the f{a) spectrum is limited to only two scales. In this sense, a two scale Cantor measure can be considered as 
a good approximation to describe the multifractal properties of natural systems. 
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